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ABSTRACT 

A grid of ATLAS9 model atmospheres has been computed, spanning 3500 K < 
Teff < 8000 K, 0.0 <\ogg< 5.0, -4.0 < [M/H] < 0.0, and -0.8 < [a/Fe] < +1.2. 
These parameters are appropriate for stars in the red giant branch, subgiant branch, 
and the lower main seq uence. The main difference from a previous, similar grid 
(ICasteUi fc Kurucj hOQl ) is the range of [a/Fe] values. A grid of synthetic spectra, 
calculated from the model atmospheres, is also presented. The fluxes are computed ev- 
ery 0.02 A from 6300 A to 9100 A. The microturbulent velocity is given by a relation to 
the surface gravity. This relation is appropriate for red giants, but not for subgiants or 
dwarfs. Therefore, caution is urged for the synthetic spectra with logs' > 3.5 or for any 
star that is not a red giant. Both the model atmosphere and synthetic spectrum grids 
are available online through VizieR. Applications of these grids include abundance anal- 
ysis for large samples of stellar spectra and constructing composite spectra for stellar 
populations. 

Subject headings: starsiatmospheres 



1. Introduction 

One approach to measuring the composi- 
tion of a star from low- or medium-resolution 
spectroscopy is to synthesize the stellar spec- 
trum. An essential ingredient in the syn- 
thesis is a model atmosphere. Typically, a 
model atmosphere tabulates quantities such 
as density, temperature, pressure, electron 
number density, and opacity clS 3j function 
of optical depth. A grid of synthetic spec- 
tra at a range of atmospheric parameters 
(effective temperature, surface gravity, and 
composition) streamlines the measurement of 
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these parameters from the observed spectrum. 
Therefore, it is useful to have a grid of model 
atmospheres, from which a grid of synthetic 
spectra may be computed. 

Mihalad (|l965l l and lStrom k AvrettI (|l965l ) 



generated the first grids of non-gray, contin- 
uum model atmospheres. Th e use of opacitv 
distr ibution functions (ODFs. lStrom &: Kurucz 
19661 ) greatly simplifies the computation of 



such grids. ODFs treat the absorption coeffi- 
cient as smoothly varying within wavelength 
subdivisions of the spectrum. In effect, the 
information about line opacity is compressed 
into a compact function, which may be used 
to compute model atmospheres. 

The grid of model atmospher es most rele- 

vant to the current work is that of Castelli &: Kurucz 
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They computed a grid of mod el atmo- 
sphere s with th e ATL AS9 program (jKurucz 
1993al ). iKurucd (|l97d ) desc ribed the original 
ATLAS code in great detail. ICastelli &: Kurucz 
sampled a wide range of effective temper- 
atures (Teff), surface gravities (logg'), and 
metallicity ([M/H]). They also sampled two 
values of [a/Fe]: 0.0 and -hO.4. 

One major use of stellar atmospheres is the 
computation of synthetic stellar spectra. Be- 
cause it is computationally expensive to gen- 
erate a model atmosphere, it is sometimes 
wise to invest the computational resources up- 
front by generating a grid of synthetic spectra. 
The red to near-infrared spectral region is the 
subject of this article. There is a precedent for 
synth etic spectral grids i n this spectral region 



(e.g.. ICoelho et al 



20051 : iMunari et al.l lioosl : 



Palacios et al.ll2010l ). The synthetic spectral 
grid here is unique because of its wide range 
in the a enhancement ([a/Fe]). 



Kirbv et al.l toid ) measured the detailed 



abundances of thousands of red giant stars 
in the dwarf satellite galaxies of the Milky 
Way. In order to do so, they required a 
grid of synthetic spectra. A grid of model 
atmospheres was computed at a range of 
[a/Fe] so that the compositions of the atmo- 
spheres would be consistent with the spectra. 
The in spiration for this proced ure was the 
grid of ICasteUi fc Kurucd (|2003l ). However, 
the presence i n dwarf galaxies o f stars with 
large f g/Fel (iFrebel et al.1 I2OI0I I and small 
[a/Fe] ( Letarte et al.ll201ol ) required a broader 
range of [ a/Fe] in the mode l atmospheres. 



Therefore, ICastelli Kuruczj 's computation 



of model atmospheres was repeated, but for 
a wide range in [a/Fe]. A grid of synthetic 
spectra was computed based on the grid of 
model atmospheres. 

The purpose of this article is to make these 
grids publicly available. Those with possi- 
ble interests in these grids include researchers 
who want to measure the elemental abun- 



dances for a large number of far-red spectra 
or who want to generate composite spectra 
of stellar populations for an integrated light 
analysis. Other uses are also encouraged. 

Table [T] gives the details of the grid param- 
eters. The step size for [M/H] is smaller for 
the spectra than for the atmospheres (see Sec- 
tion [3]). Section [2] describes the generation of 
the grid of model atmospheres, and Section [3] 
describes the grid of synthetic spectra. Sec- 
tion H] explains how the data may be retrieved 
from the VizieR online catalog. 

2. ATLAS9 Model Atmospheres 

The grid (Table [T]) of plane-parallel model 
atmospheres in local thermodynai nic equ i- 
librium (LTE) was computed with iKuruczj 's 
(Il993al l DFSYNTHE, KAP PA9, and AT - 



LAS9 codes port e d into Linux (jSbordone et al. 



20041 : ISbordonel I2OO5I I and compiled with 
the Intel Fortran compiler 11.1. All com- 
putations were performed in serial on 256 
processors at the University of California 
Santa Cruz Pleiades Linux computing clus- 
ter. The compositions of the ODFs were the 
scaled solar abundances of lAnders &: Grevesse 
(Il989l ) except that the abundance of iron was 
12 + loge(Fe) = 7.5fl The a elements O, Ne, 
Mg, Si, S, Ar, Ca, and Ti were additionally 
scaled by the parameter [a/Fe]. 

The starting point for an ATLAS9 model 
atmosphere is an ODF. O DFs were computed 
in a manner identical to Castelli &: Kurucz 



^e(X) is the ratio of the number density (n) of atoms of 
element X to the number density of hydrogen atoms: 
e(X) = n(X)/n(H). However, note that the ATLAS 
code accepts abundances expressed as n(X)/ntotai- 
This work assumes t he solar hehum-to-h ydrogen num- 
ber ratio of 0.098 () Anders fc Grevessd [TqSQ'I. Also 
note that the solar abund ance pattern assumed here 
(|Anders fc Grevessd 1 19891 ) is different from the so- 
lar abundance pattern of iGrevesse fc Sauvall (|l998h . 
which is the defauh for the DFSYNTHE and ATLAS9 
codes. 
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Table 1 
Atmospheric Parameter Grid 



Quantity Minimum Maximum Step Number 



log5 (cm Teff < 6800 0.0 5.0 0.5 11 

log 5 (cm s-2), Teff > 7000 0.5 5.0 0.5 10 

[M/H] (atmospheres) -4.0 0.0 0.5 9 

[M/H] (spectra) -4.0 0.0 0.1 41 

[a/Fe] -0.8 +1.2 0.1 21 



Note. — Sections [2] and [3] give the parameterization of microturbulent velocity, ^. 



(|2003l ) except for the elemental composi- 
tion. In summary, DFSYNTHE (described by 
Castellilliooil ^ was used to generate an ODF 



at each point in the stellar parameter grid. 
F. Castelli's scripts (http : //wwwuser . oat . ts , 



were employed for this purpose. These scripts 
compute ODFs at lower and higher effective 
temperatures than the grid presented here, 
but keeping the file formats consistent with 
Castelli's scripts simplified later steps in com- 
puting the model atmospheres. The lists of 
atomic and some molecular (including H2, 
HD, CH, C2, CN, CO, NH, OH, MgH, SiH, 
and SiO but not including TiO or H2O) tran- 
sitions were those of R. L. Kurucz, down- 



vective overshooting was turned off, and the 

mixing length parameter for convection was 

I /Hp = 1.25. For each combination of Tefj, 

log 5, [M/H], and [a/Fe], there are two at- 

astroni(lg|]dffirteliifiti) two different values of mi- 

The values are the 
1 



croturbulent velocity, ^. 
two velocities among 0, 1, 2, and 4 km s" 
that bracket the microturbulent velocity ap- 
propriate for_tlie_star'ssur^ gravity (Equa- 
tion 2 of iKirbv et aD jioO^)): ^ (km s~^) = 
2.13-0.23 log 5. 

ATLAS9 requires an input model atmo- 
sphere as an initial guess. An ATLAS9 model 
atmo sphere with similar Tpff, log 5, and [M/H] 

from Castelli Sz Kuruc^ 's grid was used for 

loadedat http : //wwwuser . oat . ts . astro . it/caste|gii/^^|$f|j(5^s/^^^lg^. ii^irions were com- 



See iKurucz. (199Q ) for details of the line lists. 
The molecular line lists for TiO and H2O, 
do wnloaded a .t the sam e URL, were provided 
bv Schwenke ( 19981 ) and Partridge Sz Schwenke 
ilMZl), respectively. The KAPPA9 code in 
the same software package as DFSYNTHE 
computed Rosseland mean opacities from the 

ODFs. 

For consistency with ICastelli &: Kuruczi 's 
(|2003l ) grid of ATLAS9 atmospheres, con- 



puted for each stellar atmosphere. After 30 it- 
erations, convergence was tested by examining 
the difference in the flux and flux derivative 
between the last two iterations. Infrequently, 
these differences exceeded 1% for the flux or 
10% for the flux derivative. In these cases, 
a different initial model atmosphere was used 
as input for 30 new ATLAS9 iterations. If the 
resulting atmosphere was still not converged, 
additional iterations were computed until the 
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7^^3500 K, logg=0.0, 
[M/H]=-1.0, [a/Fe]=0.0 



r,fp4500 K, logg=1.5, 
[M/H]=-1.0, [a/Fe]=0.0 



r,fp5500 K, logg=3.0, 
[M/H]=-2.5, [a/Fe]=0.4 




i 



O 



this work 
Castelli & Kurucz (2003) 




Fig. 1. — ATLAS9 model atmospheres computed he re (black) compare d to th e ATLAS9 model 
atmospheres with the same parameters computed by Castelli Sz Kurucz ( 20031 . gray). The black 
and gray lines overlap for much of this figure. Each column shows a different star whose parameters 
are given at the top of the column. The five rows show the mass depth (J pdx), temperature (T), 
gas pressure (-Pgas)) electron number density (ng), and Rosseland mean opacity (kr) as a function 
of the logarithm of optical depth (logr). 
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convergence criteria were satisfied. 

In very rare cases, it was not possible to 
achieve these criteria with any number of it- 
erations. These atmospheres were left as is. 
Most of these atmospheres corresponded to 
very luminous, warm red giants (6000 K < 
Teff < 8000 K and 0.0 < logg < 0.5). 
These stars lose mass in winds, which violate 
ATLASQ's assumption of hydrostatic equilib- 
rium. Stars with these atmospheric parame- 
ters are exceptionally rare. Therefore, conver- 
gence problems in this region of the grid are 
unlikely to cause problems for real astrophys- 
ical applications. The model atmospheres for 
M dwarfs (Teff < 4000 K, log 5 < 3.5) also vi- 
olated the convergence criteria to a small de- 
gree in the outer layers, especially at [Fe/H] < 
— 1. M dwarf atmospheres have low electron 
densities and non-blackbody spectral energy 
distributions, which c ause them to be out of 
statistical equilibrium ( Schweitzer 19991 ). The 
outer layers affect the strongest absorption 
lines, which are often formed out of LTE even 
for warmer and larger stars. 

Figure [1] shows some example stellar at- 
mospheres. The parameters for each atmo- 
sphere were chosen to lie near a 14 Gyr 
isochrone from the Victoria-Regina set of 



isochrones (jVandenBerg et al . 2006). For all 
three models, the microturbulent velocity is 
^ = 2 km s~^. The parameters plotted are 
those relevant to the generation of synthetic 
spectra. From top to bottom, they are the 
logarithms of the mass depth (J pdx), tem- 
perature (T), gas pressure (Pgas), electron 
number density (ue), and Rosseland mean 
opacity (kr). The independent variable is 
the logarithm of t he o ptical depth. The 
Castelh fc Kurucd (|2003l ) models for identi- 
cal stellar parameters are also shown in gray 
for comparison. Subtle differences between 
the two sets of atmospheres — most visible in 
the middle column of Figure [T] — occur due 
to the slight differences in elemental compo- 



sition and different criteria for convergence. 
The largest differences occur in the shallow- 
est layers of the atmosphere, which are impor- 
tant only for the formation of the strongest 
saturated lines. 

3. MOOG Synthetic Spectra 

The grid of synthetic spectra was computed 
with the 2007 vers ion of the pl ane-parallel, 
LTE code MOOG (ISnedenI [1973 ) . The code 
was modified to run in parallel via OpenMP 
on 256 processors on the Pleiades cluster. The 
parallelization was trivial, with different re- 
gions of the grid being computed simultane- 
ously. Each spectrum was computed on a sin- 
gle processor. MOOG was also modified to 
output only the synthetic spectrum into an 
unformatted binary file. All other output was 
suppressed to save disk space and file access 
time. The normalized flux is computed every 
0.02 A from 6300 A to 9100 A. The units of 
the spectrum are such that the continuum is 
unity. The continuum shape is not computed. 
One spectrum is computed at each grid point. 

MOOG requires two data tables as input: a 
model atmosphere and a line list. The model 
atmospheres were provided by the ATLAS9 
atmospheres described in Section [2j However, 
the spectral grid was five times finer in the 
[M/H] dimension than the model atmospheres 
(see Tabled]). Therefore, intermediate atmo- 
spheres spaced at 0.1 dex in [M/H] were gen- 
erated by linear interpolation of the physical 
variables (not the logarithms of the physical 
values) . 

The microturbulent velocity chosen for 
each synthetic spectrum is a val ue appropri 



ate f or red giants (Equation 2 of iKirbv et al 
2OO9I ): C (km s"^) = 2.13 - 0.23 logg. The 
atmosphere used to generate a spectrum was 
a linear interpolation between the two brack- 
eting values of This linear interpolation 
was performed in conjunction with linear in- 
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Fig. 2. — Examples of three synthetic spectra computed with MOOG. The stellar parameters and 
model atmospheres are identical to those shown in Figure [TJ 



terpolation in [M/H], where necessary. 
Kirbv. Guhathakurta. Sneden ( 20081 ) g; 



ave 



the complete line list, but the ve rsion used 



here h as been modified according to lKirbv et al 
(|2009l V The original line list contained 



atomic lines from the Vienna Atom ic Line 
Database fVA LD. iKuuka et al.l [l999l^. molec- 
ular lines from Kurucz J 19921). and hyperfine 
atomic lines from Kuruczj ( 1993bl ) . Oscil- 
lator strengths were modified to match ob- 
served lme__strengths_jn Arcturus and the 
Dun. iKirbv et ahf (|2009l ) replaced the oscil- 
lator strengths (I qg.q/) for soin e Fe I lines 
with the values of iFuhr fc Wiesd (120061'). For 
those Fe I lines not shared with lFuhr Wiese 
(|200(il ). 0.13 dex was subtracted from log 5/. 
The only molecules in the line list are CN, 
C2, and MgH. The opacity for almost every 
line is considered only within 1 A of the line's 
central wavelength. Some lines are considered 
at every wavelength in the spectrum. These 
are Ha, the hydrogen Paschen series redward 
of 8370 A, the Ca II triplet at 8498, 8542, and 



8662 A, and Mg I A8807. 

Figure [2] shows 100 A (7900-8000 A) of the 
synthetic spectra computed from the model 
atmospheres shown in Figure [TJ The full 
spectral range is 63 times what is shown. 
The abundances used to compute the line 
strengths are completely consistent with the 
abundances used to compute the model at- 
mospheres. 

Figure [3] shows synthetic spectra compared 
to spectra observed with the Ke ck/DEIMOS 



mediu m-resolution spectrograph. iKirby et al 



(|20ld ) obtained these spectra. These stars 
were chosen to have atmospheric parameters 
fairly close to those shown in Figures [1] and 
[2j The stars are red giants in the globular 
clusters M5 and M15. The atmospheric pa- 
rameters for the synthetic spectra are given in 
the legend at the left of each panel (measured 



bv lKirbv et al1l20ld l. The synthetic spectra 



were linearly interpolated, pixel by pixel, to 
reproduce the measured atmospheric param- 
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Fig. 3. — Examples of three synthetic spectra (black) compared to spectra observed with DEIMOS 
(gray). The synthetic spectra were interpolated in the grid in order to replicate the atmospheric 
parameters given in the legend at the left of each panel. The interpolated spectra were smoothed to 
match the resolution of the observed spectra (R ~ 6500). Some of the prominent stellar absorption 
lines are labeled with the abbreviation of the responsible element. 



eters. Finally, the synthetic spectral resolu- 
tion has been degraded to match the DEIMOS 
instrumental resolution. In the spectrograph 
configuration used, profiles of unresolved lines 
have FWHM ~ 1.2 A almost independent of 
wavelength. Therefore, the resolving power is 
R « 6200 at 7400 A and 6700 at 8700 A. 

The grid of synthetic spectra has several 
shortcomings, some of which are enumerated 
here. 



1. The only molecules included are CN, 
C2, and MgH. The number of TiO lines 
makes including them in the spectra 
computationally prohibitive. Therefore, 
cool, metal-rich stars are not modeled 
well in this grid. In reality, their spec- 
tra have TiO features that become ex- 
tremely strong as the temperature de- 
creases and metallicity increases. 



2. Strong lines are not modeled well, 
mostly due to non-LTE effects. For ex- 
ample. Figure [3] clearly shows that the 
observed wings and core of the Ca II 
A8662 line are much stronger than in 
the synthetic spectrum. 

3. Although a relation between microtur- 
bulent velocity and surface gravity has 
been established, not all stars obey 
this trend. In fa ct, the relation that 
Kirbv et al.1 (|2009l l adopted was based 
only on red giants. The relation is 
not appropriate for most subgiants and 
dwarfs. The user is strongly cautioned 
against using any synthetic spectra with 
logg > 3.5 without verifying that the 
microturbulent velocity is appropriate 
for the specific star under consideration. 
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4. Online Data Access 

The grids of model atmospheres and syn- 
thetic spectra may be obtai ned at the VizieR 
online catalogue service (lOchsenbein et al 



200d ) at the following URL: http : //vizier . u- 

The user has two options to access the cata- 
log. First, the catalog may be accessed either 
through the VizieR online interface, which in- 
cludes an option to plot the synthetic spectra 
in a web browser. Second, one or both of the 
grids may be downloaded via FTP or HTTP. 
The compressed file sizes are 310 MB for the 
atmospheres and 46 GB for the spectra. 
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